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Abstract 

We consider Lambda-hyperon production in the target-fragmentation region of semi-inclusive 
deep-inelastic scattering within the framework of fracture functions. We present a first attempt 
to determine the flavour and energy dependences of these non-perturbative distributions through 
a simultaneous QCD-based fit to available neutral- and charged-current semi-inclusive-DIS cross 
sections. Predictions based on the resulting nucleon-to-Lambda fracture functions are in good 
agreement with data and observables not included in the regression. The successful prediction 
of the Q 2 dependence of the Lambda multiplicity notably represents the first validation of the 
perturbative framework implied by fracture functions. 
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1 Introduction 



It has been known for a long time [T] that in hadronic collisions the longitudinal momentum 
spectrum of particles produced in target fragmentation crucially depends on the difference of the 
valence-parton composition of the initial- and final-state particles. In particular, only an initial- 
state particle whose valence-quark flavour content is almost or totally conserved in the scattering 
can be a leading particle in the final state, i.e. carry a substantial fraction of the incoming pro- 
jectile energy. Simple counting rules connect the differences in the valence-parton composition of 
initial- and final- state particles with the shape of longitudinal momentum spectrum of the latter. 
Leading particles in the final state are typically characterised by large longitudinal momentum 
fractions and very small transverse momenta with respect to the collision axis, a typical regime 
dominated by soft QCD dynamics where standard perturbative techniques cannot be applied. 

The leading particle effect manifests itself in reactions involving at least one hadron in the 
initial state. Quite interestingly, it also appears in processes which involve point-like probes, 
such as Semi-Inclusive Deep Inelastic Scattering (SIDIS) and, as a special case of the latter, find 
spectacular realisation in diffractive DIS. At variance with other processes discussed above, such 
process involves a large momentum transfer which enables the use of perturbative QCD. Semi- 
Inclusive DIS, in fact, provides detailed information on both the partonic structure of hadrons 
and on the fragmentation of partons into final-state hadrons. This information is generally ex- 
tracted from the so called current-fragmentation region, i.e. the phase-space region in which the 
struck parton hadronises. The target-fragmentation region is instead sensitive to the hadronisa- 
tion properties of the coloured spectator system which results from the removal of one parton 
from the incident nucleon by the virtual probe. For these reasons, the description of particles 
production in this particular region of phase space through standard perturbative calculations 
based on parton distributions and fragmentation functions will fail when compared to data. The 
description can be improved only with the introduction of new non-perturbative distributions 
which encode these peculiar aspects of QCD dynamics. 

This was early realised in Ref. [2] where the authors introduced the concept of fracture 
functions. These distributions simultaneously encode information about the parton interacting 
with the hard probe and about the fragmentation properties of the spectator system. Although 
intrinsically of non-perturbative nature, the scale dependence of such distributions can be cal- 
culated within perturbative QCD [2]. Fracture-functions evolution is driven by inhomogeneous 
Altarelli-Parisi equations which result from the structure of collinear singularities in the target- 
fragmentation region [2J [3] . Moreover, the factorisation theorem [U |5] guarantees that fracture 
functions are universal distributions, at least in the context of SIDIS, and this fact constitutes 
a solid support for phenomenological analysis. 

The phenomenology which makes use of all these concept is confined to the study of hard 
diffraction in Deep Inelastic Scattering and, to the best of our knowledge, no analysis has so 
far been attempted for particles other than proton in the final state. It is well known, however, 
that Lambda hyperon production in SIDIS on proton is mainly concentrated in the target- 
fragmentation region since Lambdas show a significant leading particle effect. 

In this paper we describe how a variety of Lambda leptoproduction cross sections can be 
simultaneously described within the fracture-function approach, if these non-perturbative dis- 
tributions are modelled at some low scale and if their free parameters are determined by a fit 
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to the available data. The rather scarce Lambda leptoproduction data in general do not allow 
to directly verify the leading-twist nature of this type of processes, which is implicitly assumed 
in the fracture-function framework, nor allow to test the scale dependence embodied by their 
specific evolution equations. In this respect, both the formalism and the model presented in this 
paper require more experimental information for a conclusive validation. We believe, however, 
that a quantitative tool that is able to reproduce many aspects of the existing data may further 
stimulate both the theoretical and the experimental activity. As a by-product, the model will 
give us the first insights on the flavour and energy dependences of the fragmentation properties 
of the spectator system into Lambda hyperons. 

The paper is organized as follows. In Section |2] we briefly recall basic SIDIS cross sections 
and fracture-function properties. In Section |3] we discuss specific features of the data sets and 
observables used in the analysis. In Section H] we describe a simple model for Lambda fracture 
functions and in Section [5] we provide and discuss the results of our fit. In Section [6] we compare 
the resulting model predictions for data and observables not used in the fit. We finally collect 
our conclusions in Section [7J 



2 Semi-Inclusive DIS in the target-fragmentation region 

The deep inelastic scattering cross section of a lepton I off a proton p with four-momenta k and 
P, respectively, is described in terms of the lepton variables: 

Q 2 P-q Q 2 n2 2 

x B = 7— — , y = -^-r = , Q =-q, (i) 

IP ■ q P ■ k ShXs 

where k' and q = k — k' are the outgoing lepton and virtual photon four-momenta, respectively, 
s h = (P + I) 2 is the centre of mass energy squared and W 2 = Shy(l — xb) + m>p is the invariant 
mass squared of the final state, with m 2 being the proton mass. The additional invariant [5] 

P ■ h _ E h 1 - cos h 
Zh ~ P-q~ E P (l-x B ) 2 ' [Z) 

is often used to specify the kinematics of final state hadron with four-momentum h, where Eh and 
Oh are the detected hadron energy and angle respectively defined in the virtual photon-proton 
centre of mass frame. The variable Zh is however not adequate to describe target fragmentation, 
since both soft hadron production (Eh — 0) and hadron production in the target-remnant 
direction (Oh — 0) both yield a vanishing value of Zh- We therefore consider cross sections 
differential either in the scaled hadron energy variable z or ( [5] 

C C=^, (3) 



1 — xb Ei 

again defined in the 7*p centre of mass frame. The final state hadron is detected in the final state 
with a fraction z G [0, 1] of the spectator energy E p (l — x B )- In the following we will analyse 
data presented in term of Feynman's variable 
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We adopt the convention that, in the 7*p frame, negative values of xf are attributed to final state 
hadrons moving parallel to incoming proton direction. Hence we assume that, in the quark-parton 
model, the target-fragmentation region is defined by the frame-dependent condition xp < 0. In 
eq. (jl]) we have introduced the hadron transverse mass, = m 2 h + p\ ± , defined in terms of its 
transverse momentum and mass squared. 

In the quark-parton model, the neutral-current semi-inclusive DIS cross section for producing 
an unpolarised Lambda off a proton in the target-fragmentation region reads [3] 



i in 



dx B dQ 2 d( Q 4 



JY + < [m£ /p (xb, Q\ C) + Mf> p (x B , Q\ 0] , (5) 

q=u,d,s 



where Y+ = 1 + (1 — y) 2 . The cross section has been re-expressed for later convenience in term 
of the ( variable, and the jacobian J = ([(1 — xb)\xf\]~ 1 has been explicitely indicated [6]. 
The latter reduces to unity in the high-energy limit and it is therefore often omitted in the 
literature. The neutrino- and anti-neutrino-induced charged-current semi-inclusive cross sections 
read respectively 

1|_ = ^J 8m , [ 2 (M^ + M*, + 2(1 - yfM^} , (6) 

and 

where the dependencies of Mi appearing in eq. (jS]) and eq. (J7J) are to be understood as in eq. ([5]) 
and the factor rjw is defined in terms of the Fermi constant Gf, the H^-boson mass and the 
electromagnetic coupling constant a em as [7] 



1 / G F M 2 W Q 2 



VW = 77 



2 y 4vr« em Q 2 + M 2 , 

As appropriate for a lowest-order calculation, we have assumed a vanishing longitudinal structure- 
function contribution in all formulas. We have further neglected charm quarks contributions. 

In eqs. (0 [6], [7j) the production of unpolarised Lambdas in the remnant direction is described 
by fracture functions M^ p (x,(,Q 2 ) I2j- These distributions express the probability to find a 
parton of flavour % with fractional momentum xb and virtuality Q 2 conditional to the detection 
of a Lambda with fraction ( of the incoming proton energy. The scale dependence of fracture 
functions is given by the following evolution equations [2] 



9 ^'(%(,fl = t 5 /' *j?'mm;"(^,( i( ,' 

2vr J Xb /(i-q u \ u 



d log n 2 

2tt J Xb x b (1-u) \ u J \x B (l-u^ 

where and P^' l (u) are the regularised [8] and real [9] Altarelli-Parisi splitting functions, 

respectively. Eq. fl9]) describes both processes which contribute to Lambda production in the 
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Table 1: Data sets used in the present analysis. {Ei) is the average energy of the incoming lepton. 
Average kinematics and production rates for the various data set, when available, are indicated. 



target-remnant direction. The homogeneous term on the right hand side of eq. ([9]) takes into ac- 
count the effects of collinear parton radiation by the struck parton i while the Lambda originates 
from the fracture function itself. The inhomogeneous one instead takes into account the possi- 
bility that observed Lambda results from the fragmentation of the radiated parton /, emitted 
collinearly to the incoming parton j. This term in fact is a convolution of parton distributions 
fj/p( x B,Q 2 ) and fragmentation functions D^(z,Q 2 ). 

For the purpose of this paper, we use leading-order formulas both for semi-inclusive cross 
sections and for the fracture-function evolution equations. We note, however, that the full formal- 
ism is available at next-to-leading-order accuracy for both unpolarised [3] and polarised [10l [TT] 
processes. 

3 Data sets and Observables 

The data used in the present analysis come from a variety of fixed-target experiments. We 
include neutrino and anti-neutrino DIS data which are crucial in providing minimal quark- 
flavour discrimination. In particular, stringent cuts (xb > 0.2 and W > 4 GeV) in data from 
Ref. [12] enhance the contributions from the valence content of the nucleon; data from Ref. [13] 
are expected to be an admixture of valence and sea contributions and therefore constrain the 
relative normalisations of the various fracture functions. We further include in the fit neutral- 
current DIS data at higher beam energy presented in Refs. [HI [15] in order to provide the 
necessary information about the energy dependence of the cross sections. 

Some care should be used in the selection of the targets. Since charged-current cross sections 
are significantly lower than neutral-current ones, many experiments have used nuclear targets. 
The first consequence is that nuclear corrections to fracture functions should be accounted for. 
Even more important for the purposes of this analysis, the fragmentation process itself might 
be affected by the nuclear medium, both in the current [16] and in the target [17] fragmentation 
region. In particular, it has been recently reported in Ref. [T7j that strange-particle yields in 
neutrino-nuclei interactions are enhanced in the target region, probably due primary particles 
re-interaction with the nuclear medium. Since the particle yields show a mild power dependence 
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on the atomic number of the target, we include in the present analysis only proton- and deuteron- 
target data. 

We would like to mention two effects which might affect the absolute normalisations of the 
various data sets. The first one is related to the definition of fracture functions. Many analyses 
quoted in Tab. ([T]) have estimated the contributions to the Lambda yield coming from the decay 
of higher-mass resonances. When not otherwise stated in the original publications, we interpret 
the published data as referring to an inclusive Lambda sample, that is the sum of promptly 
produced Lambdas and Lambdas coming from the decay of higher-mass resonances, corrected 
for unseen decay modes (a typical example is the S° — > A7 decay mode discussed in Ref. [14]). 
The decay of higher-mass resonances infact happens at time scales much larger than the ones 
typical of perturbative processes and their effects will be effectively incorporated in fracture 
functions. 

The second problem is posed by the contamination of Lambda yield by secondary Lambdas 
produced by the interaction of primary pions with detector material. This effect has been in- 
tensively studied in Ref. [T5] and estimated to contribute up to 20% to the Lambda yield. It is 
unknown to us to which extent this correction has been properly estimated and applied to all 
data sets. 

We close this Section discussing the choice and the reconstruction of the observable to be 
used in the fit. By definition, cross sections differential in the Lambda fractional energy are 
insensitive to the phase space region in which the latter has been produced. In this variable, the 
current- and target-fragmentation contributions overlap and the extraction of the latter therefore 
crucially depends on the precision with which we describe current fragmentation with available 
Lambda fragmentation functions. In order to overcome this problem, we consider differential 
cross sections in the Feynman variable xp which offer, to lowest order, a kinematical separation 
of the two contributions. The use of such a variable is however not free from additional issues: 
Lambda-mass effects introduced via eq. (0J may be sizeable, as suggested by the values of the 
averaged hadronic final-state invariant mass (W 2 ) quoted in Tab. CQ). Such effects are however 
not compatible with the pQCD factorisation theorem. In the present analysis, Lambda-mass 
effects are therefore applied a posteriori to the Lambda leptoproduction cross sections a A , as 
described in Ref. [19]. The value in each xp-bm is calculated as follows 



dE&{E{) dx B dQ 2 dC -UL &( XF ), (10) 

J Jq Jo dx B dQ 2 d( 



0"dis A-Xp 0"dis 



where the index i labels the i-th bin and the bin-size is specified by A with 
x F representing the experimental bin-edges. Mass corrections are enforced with the kinematical 
constraint l (xp) = 6{xp — x F ) 9{x l p 1 — Xp), with xp calculated via eq. (jlj). The label f2 stands 
for the set of cuts which define the DIS selection of a given analysis. The resulting differential 
cross sections are then normalised with respect to the inclusive DIS cross section 



a DIS = / dE&iEi) / dx B dQ 2 , \'l , (11) 



dx B dQ 2 

calculated with parton-distribution functions of Ref. [20] . Both cross sections are integrated over 
the lepton flux factor $(Ei) expressed in units of GeV -1 . For monochromatic electron and muon 
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beams of energy Ei t o, the latter simply reduces to 5(Ei — Ei$). For neutrino and anti-neutrino 
beams we use the flux-factor parametrisations extracted by dedicated analyses [21]. We finally 
note that mass-corrected distributions are derived by using eq. ([3]) and eq. (0| and therefore in 
principle require the knowledge of the Lambda transverse momentum. From the very precise 
data of Ref. [18J we know that the Lambda ^-spectrum is dominated by p t values much smaller 
than its mass. We therefore approximate in eq. (jl]) the transverse mass with Lambda mass 
m\ whose value is taken to be m\ = 1115.683 MeV [7]. 



4 Modelling Lambda Fracture Functions 

The fracture formalism relies on the assumption that leading particle production in the target- 
fragmentation region is a leading-twist process. The latter has been strikingly confirmed by 
experimental observation of hard diffraction at HERA [221 122] ■ in the present context, however, 
the limited amount of data is often presented as Q 2 -integrated single-differential distributions, 
a fact which prevents any conclusion to be drawn. We only note that such hypothesis might be 
indirectly supported by the moderate value of the average Q 2 for the semi-inclusive reactions 
quoted in Tab. ([T]) and by the mild increase of the Lambda average multiplicity (n(A)) as a 
function Q 2 , as seen in data [18]. 

Although the scale dependence of fracture functions is predicted by perturbative QCD, they 
still need to be modelled at some low scale and evolved to scales relevant to the experiments; the 
free parameters controlling their input distributions can then constrained by a fit to data. The 
input distributions are in general surfaces in the (xb,C) space and contain a large number of free 
parameters needed to control the energy and flavour dependences of the cross sections. Since the 
hard scattering process occurs on time scales much shorter than spectator-fragmentation ones, 
we assume that, at an arbitrary low scale Qq, fracture functions factorise into ordinary parton 

fi/p{ 



distributions fi/ p (xB,Qo) an d what we adress as spectator-fragmentation functions Df' p ' 



(1 - x B ) Mt /P {x B , £ Ql) = Mt lp {x B , z, Ql) = f i/p (x B , Ql)Df /p (z) ,i = q,q,g. (12) 

We take advantage of the sea-valence decomposition offered by the parton distributions of 
Ref. [20J to further decompose the valence parton contributions as 

K=u4^b, z, Ql) = q val (x B , Ql)D^{z) + q sea (x B , Ql)D%? a { z ) • ( 13 ) 

For the non-perturbative spectator-fragmentation function we choose a simple functional form 
of the type 

Df /p {z) = W iZ ai (1 - z)^ . (14) 



In order to minimise correlations between parameters, the normalisation coefficients in eq. ([14 
are defined as follows 



Ni = Ni 



dz z^ (1 - z) pi , cti,Pi>-l, (15) 



n -1 







and Ni are used as free parameters in the fit. 
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The inclusion of deuteron-target data in the fit requires the knowledge of neutron-to-Lambda 
fracture functions. As a first approximation we assume u-d isospin symmetry 

M^ /n (x B , z, Q 2 ) = M^(x B , z, Q 2 ), M^ n (x B , z, Q 2 ) = M^ /p (x B , z, Q 2 ) . (16) 

Limiting ourselves to the discussion of the valence region and indicating in parenthesis the flavour 
structure of the spectator system, eq. (1TB]) implies that 

U u(ud) - U d{ud) aIia U d{uu) - U u{dd)- 

The validity of these assumptions will be discussed in Section |6j 

Given the input distributions in eqs. ( TT21T3I) . set at an initial scale Ql = 0.5 GeV^ 2 , we numeri- 
cally solve the fracture-functions evolution equations in eq. (J9j by using a finite-difference method 
in x^-space in slices of (. We use the leading-order proton parton distribution of Ref. [20] and A 
fragmentation functions of Refs. [19] extracted taking into account target-hadron mass effects. 
For consitency, we follow the original evolution scheme of Ref. [20]. We evolve light-quarks frac- 
ture functions to leading-logarithmic accuracy by using a fixed flavour-number scheme. Heavy- 
quarks effects are included in the running of the strong coupling at the heavy-quark mass thresh- 
olds and Aq^d values are taken from Ref. [2D] . 

As a final remark to this section we would like to discuss the impact of the inhomogeneous 
term in eq. (J9j on the observables entering the fit. As a first step we have checked that, in 
the current region where the leading order semi-inclusive cross sections are proportional to the 
product of parton distributions and fragmentation functions, a reasonable description of fip data 
of Ref. [U] can be obtained both in shape and in normalisation by using parton distributions 
and fragmentation functions of Refs. [20] and Ref. [19] , respectively. We then use fragmentation 
functions of Ref. [19] in the calculation of the inhomogenous term appering in eq. (Q. At the 
observable level such term accumulates, as expected, at small | result of the z shape of 

Lambda fragmentation functions. Moreover it contributes only at the percent level to the xp dis- 
tributions since it describes Lambda production from the fragmentation of radiated partons with 
transverse momentum greater than the minimum momentum transfer involved in the process, 
that is Q\. As such it only contributes to the small radiative tail of the Lambda pt- spectrum. 

The radiative contributions generated by the inhomogeneous term and the universality issue 
of fragmentation functions could be studied in detail in processes that do not show any leading 
particle effect, for example the production of anti-Lambdas [JJ] or light mesons [23J off protons, or 
in target-like Lambda production in a perturbative regime (i.e. with sufficiently large transverse 
momentum, p t > 1 GeV). 

5 Fitting procedure 

Each of the assumed Lambda spectator-fragmentation functions in eq. fll4p contains three free 
parameters (one normalisation and two exponents). Our fitting strategy is the following: we first 
tentatively assume as many different fragmentation functions as possible types of struck partons 
(valence or sea u, valence or sea d, other sea quarks, gluon). The observable l/a ms dcr A /dxF is 
then reconstructed via eq. (fTDl and the best fit parameter values are determined using MINUIT [25] 
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Table 2: Partial x 2 contributions and number of points in the fit for each data set. 
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0.046 ± 0.006 
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dval 


0.027 ± 0.006 




1.28 ± 0.51 


Qsea 


0.078 ± 0.010 





1.84 ± 0.63 



Table 3: Best-fit values according to eq. (fTil) . 



program. We assume that the uncertainties on the cross sections combine statistical and system- 
atic errors so that we use the simplest version of the x 2 -f unc ti° n as a merit function, although 
this neglects correlations between data points (the number of the latter is rather limited and 
amounts to 46). We then study the eigenvalues and the eigenvectors of the fit covariance matrix 
to identify any parameters that are badly constrained by the fit. We fix them by making some 
assumptions about their values and we repeat the fit with the remaining parameters. 

With the available data, for example, the fit can not constrain distinct sea-quark fracture 
functions. We therefore assume a common spectator-fragmentation function for all of them, D^/f a . 
Moreover the fit is found to be insensitive to the choice of the gluon spectator fragmentation 
function. This is not unexpected since the gluon is coupled to electroweak probes only through 
higher orders. For this reason, we fix the gluon spectator-fragmentation functions to be equal 
to the sea one, Dg = Dq/f a , reducing the number of free parameters to nine. In such a fit the 
smallest eigenvalues of the Hessian matrix correspond to eigenvectors whose largest components 
are associated with the parameters «j. The poor determination of such parameters can be partly 
associated to mass effects in the reconstruction of the observable in eq. (|T0|) via eq. (@J. The 
parameter a 0gea is compatible with zero within errors so we fix it to this value. Furthermore, the 
parameters a Uval and a.d val are equal to each other within errors, so that we assume a Uval = ocd val , 
reducing the number of free parameters to seven. Fixing some of the parameters to a definite value 
is indeed an arbitrary procedure. We can motivate this choice only a posteriori by noting that, 
when these parameters are optimsed by the fit, the value of the x 2 function is only marginally 
reduced. The best seven-parameter fit yields a \ 2 jd.o.f. = 44.14/(46 — 7) = 1.13. The partial 
X 2 and the number of points included in the fit for each data set are displayed in Tab. (T5]). 
The results for the best-fit parameters are reported in Tab. ([3]) along with parameter errors as 
calculated by the MINUIT routine HESSE, which assumes a parabolic behaviour of the \ 2 function 
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in parameter space around the minimum. 

The best-determined parameters are the three normalisations, which implies that three in- 
dependent distributions are in fact sufficient to handle the normalisation spread between the 
various data sets. The four parameters controlling the shape of the spectator functions have 
substantially larger uncertainties. The sources of the latter are primarily related to the intrinsic 
correlation between the on and fii parameters introduced by the specific functional form as- 
sumed for the spectator functions in eq. (TH|) . We also note that, at low centre-of-mass energy, 
where mass corrections introduced via eq. (j3J) are sizeable, a large portion of the xf spectrum 
is controlled by the behaviour of the spectator functions in the neighbourhood of z = 1; the 
shape of the predicted spectrum is thus mainly determined by the fa parameter alone. There- 
fore, parametrisations with more modulabie behaviour at large z could improve the description 
of the xp spectra. These improvements are probably marginal given the quality of the data used 
in this fit, but they might be necessary when dealing with higher quality data as presented, for 
example, in Ref. |18j . 

Given the stringent cut on the DIS selection, the data of Ref. [12J constrain the valence- 
quark fracture- function contributions, M qval , which in fact almost saturate the spectrum, as 
shown in the left column of Fig. ([T]). The first row of Fig. ([T]) reveals instead a normalisation 
tension between up data from Ref. [13] and Ref. [12] which however can be tolerated in view of 
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Figure 2: Best fit predictions compared to normalised xp distributions for /ip (left panel) and [iD<i 
(middle panel) data from Ref. [T3] and 11D2 (right panel) data from Ref. [TS]. In the \i"p case valence- 
and sea-quark fracture functions contributions are separately shown. In the 11D2 case the proton and 
neutron contributions are separately shown. 



the partial x 2 presented in Tab. ([3]). The second row shows a slight shape deformation in the 
predictions for un data from Ref. [12] which is probably due to the normalisation constraint 
induced by up data from Ref. [T3] on the individual M Uval and M qsea distributions. The fit indeed 
indicates, as expected from the valence-parton counting rules [TJ, that Lambdas are produced 
more abundantly and more forward by the fragmentation of w<i-spectator system with respect 
to a uu- one. 

We show in Fig. (j2J) the best-fit predictions for data from Ref. [H] and Ref. [15] for which the 
incident muon beam energy is significantly higher than neutrino and anti-neutrino ones. Quite 
interestingly, the spectrum at large \xp\ is dominated by valence-quark fracture functions, as 
shown in the upper left panel of Fig. ([3]). The differences between M Uval and Ma val distributions 
is responsible for the different large \xp\ behaviour on different targets, as shown in the upper 
right panel of Fig. ([3]). In all these plots, given the higher values of (W 2 ) involved, Lambda- 
mass corrections play a less prominent role and therefore the z shape of the sea-quark fracture 
functions, as parametrised in Tab. ([3]), are clearly visible in these distributions. The results 
discussed above are stable against the variation of the arbitrary scale Q\ at which fracture 
functions are factorised into parton distributions and spectator-fragmentation functions. The \ 2 
function, in fact, shows a very mild dependence on Qq as the latter is varied below the measured 
range between 0.5 and 1.0 GeV 2 . 
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Figure 3: Best-fit predictions compared to normalised Xi?-distributions for charged current semi- 
inclusive Lambda cross-sections from Refs. [26} I27j. Various quark-flavour proton-to-Lambda fracture 
functions contributions are shown. 



6 Predictions 

In this last section we wish to discuss the degree of predictivity of the model. It is of particular 
importance to determine to which extent the model is able to reproduce xf distributions coming 
from data not included in the fit or observables other than the ones used in the fit. At the same 
time, the quality of the agreement with the data will pinpoint sectors of the model which may 
need to be improved. 

The comparisons often require the reconstruction of the current fragmentation term, which is 
obtained with the semi-inclusive version of eqs. ( 161171) where fracture functions are substituted by 
appropriate products of parton distributions and fragmentation functions. Such term has been 
estimated to leading order by using the parton distributions of Ref. [20j and the fragmentation 
functions of Ref. [19] . 

As a first example, we consider Lambda production in neutrino- and anti-neutrino-induced 
charged- current DIS on proton and neutron targets [26]. These data provide a full flavour dis- 
crimination of the spectator system and can be used to test the assumptions made in eq. ( TT6l) . 
Unfortunately, the xp distributions for these data are presented in Ref. [26] as histograms with- 
out errors. In a subsequent paper by the same collaboration [27], the Lambda yields are updated 
but no xf distributions are shown. 
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Figure 4: Normalised xb distributions in vD<i (left panel) and vn (right panel) scattering. Data are 
taken from Ref. |27j . 



In order to gauge the agreement of the best-fit model predictions with these data, we have 
assessed the errors in the following way. First the histograms from Ref. [26] have been scaled 
down to the updated yields of Ref. [27] assuming no change in their shape. Then we have assumed 
that the relative errors of the xp distributions in each xp bin are equal. This in turn implies 
that they are respectively equal to the relative error on the corresponding Lambda yield. This 
procedure guarantees that upon integration over xp, the experimental yield and its error are 
correctly recovered. 

The best-fit predictions are compared to data in Fig. ([3]). We note that, in general, the cross 
sections on proton target are fairly reproduced in normalisation, as expected, since the Lambda 
yields in vp- and i^p-scattering of Ref. [27] are in agreement within errors with those of data 
included in the fit [131 [12]. The distributions show however some stress in shape at low values of 
\xp\. The shapes of the model predictions reflect those of the data used in the fit [Ml [12], which 
are shown in the right- hand- side panels of Fig. ([T]). The comparison with neutron-target data 
reveals that, in this case, the model does not perform equally well. In particular, the vn cross 
sections seem to peak at too large values of \xf\, a behaviour mainly driven by the vn data of 
Ref. [12] and illustrated in the lower left panel of Fig. While all these distributions are quite 
well reproduced in normalisation, the model significantly underestimates the z/n-scattering data. 
The level of disagreement of data and theory might indicate that D^fts < D^j^y contrary to 
our assumptions in eq. ([TBI) . This interpretation is however subject to caution, as we discuss in 
the following. 

The same data are presented in Fig. (j3J) in terms of normalised xb distributions. We remind 
the reader that we integrate over the xb variable in the reconstruction of the xp distributions 
(eq. (TTU|) ). Therefore the comparison of the xb distributions with the data constitutes a non- 
trivial test for the model. Excellent agreement is found for the vD 2 cross section, both in shape 
and in normalisation, as already seen for the vp and vp xp distributions from the same data 
set and shown in the left-hand-side panels of Fig. ([3]). On the other hand, model predictions 
for the vn cross section fail to reproduce experimental distributions, as in the corresponding xp 
distribution in the lower right panel of Fig. ([3]). Since the largest differences in the vD 2 and vn 
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Figure 5: Lambda average multiplicity as a function of Q 2 in up (left panel) and up (right panel) 
scattering. Current and target contributions are separately shown. Data are from Ref. |13j . 



cross sections appear in the smallest xb bin, there might exist some asymmetry in the sea-quark 
spectator functions which we instead assume to be all equal to each other. 

It was originally suggested in Ref. [26J that target Lambdas are more abundantly produced 
by z/s-scattering than in i/s-scattering, based on the idea that ^-scattering leaves the correct 
hypercharge in the target spectator. Such a process would indeed enhance the yield at small xb- In 
our framework, this hypothesis should be translated into an asymmetry of the spectator functions 
in the strange sector, d£' p > Dg. While a general asymmetry in the sea-quark spectator 
functions can be easily accomodated in the fit by increasing the number of free parameters, the 
cross sections used in the regression cannot separately constrain the parameters of d£' p and 
D^ p , because of the linear combinations of sea-quark fracture functions that enter the structure 
functions in eqs. fl6|7|) . 

In this respect, leptoproduction of target Lambdas in dimuon (anti-)neutrino-nucleon deep- 
inelastic scattering could represent an alternative sensitive test of the latter hypothesis. In this 
process, a (anti-) charm, produced by a charged current, semileptonically decays in the current- 
fragmentation region into a final-state, secondary, muon. Given the small off-diagonal quark- 
mixing CKM matrix elements, this process is able to probe directly the (anti-)strange component 
of the nucleon [28] and consequently the proposed Lambda asymmetry, > M^ p , in the 

target-fragmentation region. 

To conclude we show in Fig. (J5J) the model predictions for the averaged Lambda multiplicity 
as a function of Q 2 in up and up charged-current DIS cross sections. Although the xp distri- 
butions from this data set [13] have been already used in the fit and therefore agreement in 
normalisation might be expected, the model is also able to reproduce reasonably well the ob- 
served Q 2 dependence, with a tendency to undershoot the data at the lowest value of Q 2 . Since 
the Q 2 dependence built in the model via fracture-function evolution equations can be considered 
one of the most stringent predictions of the underlying theoretical framework, the reasonable 
agreement between data and predictions can be considered the first step towards a conclusive 
validation of the perturbative framework. The validation procedure and the model itself would, 
in fact, highly benefit from the constraints imposed by multi-differential distrubutions. Cross 
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sections at fixed xb and function of Q 2 , for example, could give access, through scaling 

violations, to the presently unconstrained gluon fracture function, M^ p (xb, z,Q 2 ). 

7 Summary and conclusions 

In this paper we have analysed experimental data on the production of Lambda hyperons in 
the SIDIS target-fragmentation region in terms of fracture functions. A model for the latter has 
been proposed and the free parameters appearing in the input distributions have been fixed by 
performing a fit to a variety of neutral- and charged-current semi- inclusive DIS cross sections. 
The main features seen in the data can be fairly reproduced by the model. In particular the 
spectator-fragmentation functions associated with the removal of valence quarks populate the 
very forward part of the Xi?-spectrum at large and negative values of xp. On the other hand, the 
sea-quark contribution is concentrated at small and negative values of Xp. The predictions based 
on the model are in fair agreement with data not included in the fit and with observables depend- 
ing on variables which are integrated over in the analysis, especially the Q 2 dependence of the 
Lambda multiplicity, which is a stringent test of the underlying theoretical framework. Although 
further tests and additional experimental informations are necessary to validate and eventually 
improve the model, it may be used to quantitatively investigate spectator-fragmentation mecha- 
nisms within a perturbative QCD approach. Since higher-order calculations are available in the 
literature, the analysis can be extended to next-to-leading order accuracy. The model can be 
easily generalised to take into account Lambda polarisation allowing spin-transfer studies [29] in 
the target region and it may find application in the estimation of nuclear corrections to target 
fragmentation. 
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